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Abstract. First we introduce and analyze a convergent numerical method 
for a large class of nonlinear nonlocal possibly degenerate convection diffusion 
equations. Secondly we develop a new Kuznetsov type theory and obtain 
general and possibly optimal error estimates for our numerical methods — even 
when the principal derivatives have any fractional order between 1 and 2! The 
class of equations we consider includes equations with nonlinear and possibly 
degenerate fractional or general Levy diffusion. Special cases are conservation 
laws, fractional conservation laws, certain fractional porous medium equations, 
and new strongly degenerate equations. 



1. Introduction 

In this paper we develop a numerical method along with a general Kuznetsov 
type theory of error estimates for integro partial differential equations of the form 

f d t u + div/(«) = [A(u)} , {x, t) e Qt, 
( ' \ u(x,0) = u (x), x e M d , 

where Qt = M. d x (0,T) and the nonlocal diffusion operator £ M is defined as 



(1.2) C[4>]{x)= 4>(x + z)-4>(x)-z>V<Kx)l\ t \ <1 (z)dn(z), 

J\z\>Q 

for smooth bounded functions <p. Here 1 denotes the indicator function. Through- 
out the paper the data (/, A, (j,,uq) is assumed to satisfy: 

(A.l) / = {f x , ...J d )£ W 1 ' 00 (R;R«') with /(0) = 0, 

(A.2) A e W 1 '°°{M), A non-decreasing with A(0) = 0, 

(A. 3) ii > is a Radon measure such that J| z | >0 \z\ 2 A 1 dfi{z) < oo, 

(A.4) u G L°°(R d ) n L^W 1 ) n BV{R d ). 
We use the notation aAi) = min(a, b) and a V b = max(a, b). 

Remark 1.1. These assumptions can be relaxed in two standard ways: (i) /, A can 
take any value at u = (replace / by / — /(0) etc.), and (ii) /, A can be assumed 
to be locally Lipschitz. By the maximum principle and (A.4), solutions of (|1.1[) are 
bounded, and locally Lipschitz functions are Lipschitz on compact domains. 

The measure \i and the operator D 1 are respectively the Levy measure and the 
generator of a pure jump Levy process. Any such process has a Levy measure 
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and generator satisfying (|1.2[) and (A. 3), see e.g. [I]. Example are the symmetric 
a-stable processes with fractional Laplace generators where 
dz 

(1.3) dii(z) = c xrm z (c A >0) and = -(-A) A / 2 for A e (0,2). 

PI 

Non-symmetric examples are popular in mathematical finance, e.g. the CGMY 
model where 

C e -G| 2 | 

— — — -rr — dz for z > 0, 

(j e -M\z\ 

——r—rr—dz for z < 0, 

|z| i+A 

and where d = 1, A(= Y) G (0, 2), and C, G, M > 0. We refer the reader to [H] for 
more details on this and other nonlocal models in finance. In both examples the 
nonlocal operator behaves like a fractional derivative of order between and 2. 

Equation (11 . If) has a local non-linear convection term (the /-term) and a frac- 
tional (or nonlocal) non-linear possibly degenerate diffusion term (the A-tevm). 
Special cases are scalar conservation laws {A = 0), fractional and Levy conserva- 
tion laws (A(u) = u and a-stable or more general //) - see e.g. [51 Q] and [71 1551 125] . 
fractional porous medium equations [TB] (A = \u\ m ^ 1 u for m > 1 and a-stable //), 
and strongly degenerate equations where A vanishes on a set of positive measure. 
If either A is degenerate or £ M is a fractional derivative of order less than 1, then 
solutions of (jl.ip are not smooth in general and uniqueness fails for weak (distri- 
butional) solutions. Uniqueness can be regained by imposing additional entropy 
conditions in a similar way to what is done for conservation laws. The Kruzkov 
entropy solution theory of scalar conservation laws [27J was extended to cover frac- 
tional conservation laws in [T] , to more general Levy conservation laws in [35] , and 
then finally to setting of this paper, equations with non-linear fractional diffusion 
and general Levy measures in [llj . For local 2nd order degenerate convection dif- 
fusion equations like 

(1.4) <9 t u + div/(u) = AA(u), 

there is an entropy solution theory due to Carrillo [5]. 

In recent years, integro partial differential equations like (jl.ip have been at the 
center of a very active field of research. A thorough description of the mathematical 
background for such equations, relevant bibliography, and applications to several 
disciplines of interest can be found in [TJ [51 [7J EES HS1 125] ■ 

The first contribution of this paper is to introduce a numerical method for equa- 
tion (jl.ip and prove that it converges toward the entropy solution of (jl.ip under 
assumptions (A.1)-(A.4). The numerical method is based upon a monotone fi- 
nite volume discretization of an approximate equation with truncated and hence 
bounded Levy measure. Essentially it is an extension of the method in [TT] from 
symmetric a-stable to general Levy measures, but since non-symmetric measures 
are allowed, the discretization becomes more complicated here. Apart from its abil- 
ity to capture the correct solution for the whole family of equations of the form 
(jl.ip , the main advantage of our numerical method is that it allows for a complete 
error analysis through the new framework for error estimates that we develop in 
the second part of the paper. 

The second, and probably most important contribution of the paper, is the devel- 
opment of a theory capable of producing error estimates for degenerate equations 
of order greater than 1. This theory is based on a non-trivial extension of the 
Kuznetsov theory for scalar conservation laws [35] to the current fractional diffu- 
sion setting. An initial step in this analysis was performed in [2] , with the derivation 
of a so-called Kuznetsov lemma in a relevant form for (jl.ip . In [5] the lemma is 
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used in the derivation of continuous dependence estimates and error estimates for 
vanishing viscosity type of approximations of ([l.ip . In the present paper, we show 
how it can be used in solving the more difficult problem of finding error estimates 
for numerical methods for (|l.ip . 

As a corollary of our Kuznctsov type theory, we obtain explicit A-dependent 
error estimates when \x is a measure satisfying 

dz 

(1.5) < l\ z \ <x dn{z) < c x -rid+x for c A > and A e (0,2). 

\ z \ 

In this paper we will call such measures fractional measures. For example for the 
implicit version of our numerical method (I3.5p . we prove in Section [S] that 

( Axi Ae (0,1), 

\\u(;T)-u Ax (;T)\\ Lim <C T I Ao:5log(A a; ) A = l, 

{ Ax^ Ae (1,2), 

where u is the entropy solution of (jl.ll) and ua x is the solution of p.5[) . Note that 
our error estimate covers all values A e (0, 2), all spacial dimensions d, and possibly 
strongly degenerate equations! Also note that under our assumptions, the solution 
u possibly only have BV regularity in space. Hence the error estimate is robust in 
the sense that it holds also for discontinuous solutions, and moreover, the classical 
result of Kuznetsov [25J for conservation laws follows as a corollary by taking A = 
(a valid choice here!) and A e (0, 1). The above estimate is also consistent with 
error estimates for the vanishing A-fractional viscosity method, 

d t u + divf(u) = -Ax(-A) x/2 u as Ax -s- + , 

see e.g. [Hill], but note that our problem is different and much more difficult. 

There is a vast literature on approximation schemes and error estimates for scalar 
conservations laws, we refer e.g. to the books [55] and references therein for 
more details. For local degenerate convection-diffusion equations like (|1.4[) . some 
approximation methods and error estimates can be found e.g. in [301 HD HI] and 
references therein. In this setting it is very difficult to obtain error estimates for nu- 
merical methods, and the only result we are aware of is a very recent one by Karlsen 
et al. [24 (but see also 10 J. This very nice result applies to rather general equa- 
tions of the form (|1 .4[) but in one space dimension and under additional regularity 
assumptions (e.g. d x (A(u)) e BV). When it comes to nonlocal convection-diffusion 
equations, the literature is very recent and not yet very extensive. The paper [15) 
introduce finite volume schemes for radiation hydrodynamics equations, a model 
where £^ is a nonlocal derivative of order 0. Then fractional conservation laws are 
discretized in [T7] [T3l [H] with finite difference, discontinuous Galerkin, and spectral 
vanishing viscosity methods respectively. In [151 113] Kuznetsov type error estimates 
are given, but only for integrable Levy measures or measures like (|1.3p with A < 1. 
Both of these results can be obtained through the framework of this paper. In 
[T2] error estimates are given for all A but with completely different methods. The 
general degenerate non-linear case is discretized in [TT] (without error estimates) 
for symmetric a-stable Levy measures and then in the most general case in the 
present paper. 

Linear non-degenerate versions of (jl.ip frequently arise in Finance, and the prob- 
lem of solving these equations numerically has generated a lot of activity over the 
last decade. An introduction and overview of this activity can be found in the book 
[14] . including numerical schemes based on truncation of the Levy measure. We 
also mention the literature on fractional and nonlocal fully non-linear equations 
like e.g. the Bellman equation of optimal control theory. Such equations have been 
intensively studied over the last decade using viscosity solution methods, including 
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initial results on numerical methods and error analysis. We refer e.g. [SJ[5J[23] and 
references therein for an overview and the most general results in that direction. 
In fact, ideas from that held has been essential in the development of the entropy 
solution theory of equations like (jl.ip , and the construction of monotone numerical 
methods of this paper parallels the one in [5]. However the structure of the two 
classes of equations along with their mathematical and numerical analysis are very 
different. 

This paper is organized as follows. In Section[2]we recall the entropy formulation 
and well-posedness results for (jl.ip of [11] and the Kuznetsov type lemma derived 
in [2] . We present the numerical method in Section [3] There we focus on the 
case of no convection (/ = 0) to simplify the exposition and focus on new ideas. 
In Section 0] we prove several auxiliary properties of the numerical method which 
will be useful in the following sections. We establish existence, uniqueness, and a 
priori estimates for the solutions of the numerical method in Section [5j The general 
Kuznetsov type theory for deriving error estimates is presented in Section [51 where 
it is also used to establish a rate of convergence for equations with fractional Levy 
measures, i.e. (|1.5p holds. In Section [7] we extend all the results considered so far 
to general convection-diffusion equations of the form (jl.ip with / ^ 0. Finally, we 
give the proof of the main error estimate Theorem 16.11 in Section [8] 



In this section we briefly recall the entropy formulation for equations of the form 
(jl.ip introduced in jllj , and the new Kuznetsov type of lemma established in [2] . 
Let r)(u,k) = \u - k\, r/'(u,k) = sgn (u - k), qi(u, k) = r)'(u,k) (fi(u) - fi(k)) for 
I = 1, . . . , d, and write the nonlocal operator £ M [</>] as 



We also define fi* by H*(B) — fi(—B) for all Borel sets B ^ 0. Let us recall that 



for all smooth L°° n L 1 functions ip, ip, cf. [21 [TT]. 

Definition 2.1. (Entropy solutions) A function u g L°°{Q T ) n C([0, T];L 1 (R d )) 
is an entropy solution of (jl.ip if for all k G M, r > 0, and test functions < ip G 
C™(R d x [0, T]), 



2. Preliminaries 



where 






(2.1) 



+ rf(u,k)£^ r [A(u)]ip dxdt 




Note that 7, = when the Levy measure /1 is symmetric, i.e. when /1* = fi. 
From [TT] we now have the following well-posedness result. 
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Theorem 2.1. (Well-posedness) Assume (A.l) - (A. 4) hold. Then there exists a 
unique entropy solution u of such that 

u e L°°{Qt) n C{[Q,T];L 1 {m. d )) n L°°(0,T;BV(R d )), 

and the following a priori estimates hold 

\\u(-,t)\\L=°(Wl d ) < ||wo||l°°(R<*), 
ll tt (*)*)IU 1 (K d )) ^ ll M o|Ui(Rd), 
l u ('i*)lsU(R d ) < l u o|_BU(R<i), 
\\u(;t) - «(-,s)||ii(B<J) < — 

/or t,s£ [0, T] where 

cr I\z\>o \ z \ A 1 d ^( z ) < 00 ' 



>c 

cr' otherwise. 



a(r) 

Moreover, if also (|1.5[) holds, then 

(cr if \e (0,1), 

cr(r) = < c |r In r| if X = 1, 

[cri if Xe (1,2). 

The last a priori estimate is slightly more general then the one in |llj , and follows 
e.g. in the limit from the estimates in Lemmas 15.31 and 15.41 We now recall the new 
Kuznetsov type of lemma established in [2] . Let 

cjeC c °°(K), 0<w<l, w(r) = for all |r| > 1, and /w(r)dT = l, 

and define u>s(t) = 4 u; (?) , 17 e (x) = w e (xi) • • • w e (xd), and 

^ e,,5 (z,2/,M) = O e (x - w 5 (t - s) 

for e, 5 > 0. We also need 

£«(«)= sup \\v(-,t)-v{-,s)\\ L i {Rd) . 
(2.2) \t-s\<5 

i,se[0,T] 

In the following we let dw — dx dt dy ds and C't > be a constant depending on 
time and the initial data uq that may change from line to line. 

Lemma 2.2. (Kuznetsov type of lemma) Assume (A.l) - (A. 4) hold. Let u be the 
entropy solution of (|1.1[) and v be any function in L°°(Qt) H C([0, T]; L 1 (M. d )) n 
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L°°(0,T;BV{R d )) with v(- ,0) = v (-). Then, for any e,r > and < 5 < T. 
IK-,T) - v(-,T)\\ L i^ < \\u - uoIIl^r") +C(e + f 5 (M) V£ 5 (w)) 

?y(i>(a;, t), it(y, s)) d t tp e ' s (x, y, t, s) dw 



Qt JJQt 

q(v(x, t), u(y, s)) ■ V x cp e ' S (x, y, t, s) dw 



Qi 

[[ r,(A(v(x,t)),A(u(y,8)))Cf fe> e,a (av,M)](ff) dw 

Qt JJQt 

r,'(v(x, t),u(y, s)) C^ r [A(v(-,t))}(x) ^ s (x, y, t, s) dw 

Qt JJQt 

v(A(v(x, t)),A(u(y, a))) 7 M *' r • V x <f e > s (x, y, t, s) dw 

Qt JJQt 

/ v{v(x,T),u(y,s))(p e ' S (x,T,y,s) dxdyds 

Qt J'K. d 

Si 



r){v (x),u(y,s))ip et (x,0,y,s) dxdyds 

The proof is given in h2] . The original result of result of Kuznetsov in [28] is a 
special case when n = (or A = 0). 

3. The numerical method 

In this section we derive our numerical method. Here and in the following sec- 
tions we focus on the case / = to simplify the exposition and focus on the new 
ideas. The general case / 7^ will then be treated at the end, in Section [7] 

We will consider uniform space/time grids given by x a = a Ax for a £ Z d and 
t n = n At for n = 0, . . . , N = We also use the following rectangular subdivisions 
of space 

R a = Xa + Ax (0, l) d for aeZ d . 
We start by discretizing the nonlocal operator, replacing the measure /x by the 
bounded truncated measure l| 2 ,| > A 2 .(^)/i and the gradient by a numerical gradient 

(3.1) D Ax = (£>!,••• ,D d ), 

where Di = D? are upwind finite difference operators defined by 



(3.2) 



D+#x) := 0(X + A ;; ) '^ ) for tf* > 0, 

D7 <b(x) :— < t-^l i-i x otherwise. 

1 Aa 



\x 

Here ei, . . . , is the standard basis of M. d . This gives an approximate nonlocal 
operator 

mA(4>)](x) 

Ll4i =/ A(cj>(x + z))-A(4>(x))dfi(z)+ 7 ^ ■ D Ax A(cj>(x)), 
J\z\>4± 

which is monotone by upwinding and non-singular since the truncated measure is 
bounded. 

A semidiscrete approximation of (jl.ip with / = is then obtained by solving 
the approximate equation 

(3.4) d t u = &[A(u)], 
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by a finite volume method on the spacial subdivision {R a } a - he. for each t, we 
look for piecewise constant approximate solution 

U(x,t) = U p (t)l Rf ,(x), 



that satisfy Q3.4p in weak form with -^-^Ir^ as test functions: For every a £ Z d , 
1 /' 1 



f d t Udx = -^ f fr[A(U)]dx. 

JR a ^ x JR a 



Ax d 

Finally we discretize in time by replacing dt by backward or forward differences 
D^ t and U a (t) by a piecewise constant approximation [/". The result is the implicit 
method 



(3.5) U2 +L =V2 + AtC' l (A(U n+L )) a , 
and the explicit method 

(3.6) L/™ +1 = U£ + Atfr(A{U n )) a 
where 

t»{A{U»)) a = f fr[A(U n )](x) dx, 

and U n {x) — X^ez d Up ^-Rf>( x ) ^ s a piecewise constant x-interpolation of U. As 
initial condition for both methods we take 

U° = i ! uo(x) dx for all a € Z d . 



Lemma 3.1. 



/3GZ 

wii/i G| = G Q:(3 + G"" 3 and 

G a ,/3 = -7—7 / / 1r p (z + 2) - liis (x) d/x(z) dx, 

Ax d jR c ,J\z\>^ 

(3-7) „ ^ 

G*J> Djl R ,(x)dx. 

Remark 3.2. G Qi/ 3 is a Toeplitz matrix (cf. Lemma [4.11 (b)) while G a '^ is a tridi- 
agonal matrix. When the measure /1 is symmetric, then G a ^ is symmetric and 
G a ' = 0. 

Proof. Since 

A(tf(x)) = ^ A^l^x) and D Ax U(x) = ^ Upb/^lR^x), 
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we find that 

Ax d &{A{U)) a = / &[A{U)](x) dx 



A(U(x + z)) - A(U(x)) dn(z) dx 

R a J\z\>^f 



d \JR a J\; 



f A{Up)D Ax l Rfi {x) dx 

l Rfj (x + z)- 1r p (x) dfi(z) dx^j 
£ A(U P ) I D]l Rf! (x)dx 

3GZ d \l=l jRa / 



The proof is complete. □ 

4. Properties of the numerical method 

In this section we show that the numerical methods are conservative, monotone 
and consistent in the sense that certain cell entropy inequalities are satisfied. We 
start by a technical lemma summarizing the properties of the weights Gp defined 
in (TO) . 

Lemma 4.1. 

(*) G% = £ Q£Zd Gi = for all /3 e Z d . 

(b) G£ = G%\% for all a, /3 £ Z d and I = 1, . . . , d. 

(c) Gp <0 and G^>0fora^ j3. 

(d) There is c — c(d 7 /i) > such that G^ > — ~ fstz) an< ^ w ^ ere 

s when J \z\ A 1 d/i(z) < oo, 
s 2 otherwise. 



(4-1) a^s) 



(e) If p. 51) holds, then there is c = c(d, A) > such that Gp > — ^(Az) f or 

!s x for A > 1, 
s for A < 1. 
Proof, (a) By the definitions of G a ,p, G a '^ and Fubini's theorem, 

Ax d V G a ,p = ( l R Jx + z) dx - / l R Jx) dx] d^{z) = 0, 

Ax d V G^ = ±V^— - ( / l % (x±Aare,)d!r- / 1*, (x)dx) = 0, 
agzd l=1 Ax \J Rd J Rd J 

and, since J^pez* ^(z) = 1 > 



Ax d J2 G *,P = [ f J^Wx + z)-^ l RfJ (x)) dfi(z) dx = 0, 

^ eZ d •'fl a .'|z|>% 2 p eZ d p eId d 

Ax d £ ^^iE 2 ^ / (£Wz±AzeO-£ l^(x))dx = 0. 
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Therefore £ aeZ « G^ = J2 a ez d + G^) = and G% = 0. 

(b) Let y = x + ei and note that 

Ax d G a ^=j / l R/3 (y - Axei + z) - l Rl) (y - Axei) dfi(z)dy 
jR a+<!l J\z\>^ 

= / (v + z ) - (y) d M0) d y 

In a similar fashion we get G"'' 3 = G^ +ei,a+ei . 

(c) Note that 

Ax d Gfi.fi = f f l R Jx + z)-l d/x(z) dx < 0. 

jRp J\z\>4± 



while by the definition Dj , see (|3.2I) 



z=i 

For f3, 



Ax d G^ = -V 7 r ¥ sgn( 7 f'^) / l^M da; <0. 



Ax d G a .fi= / l R (x + z) d/i(z)dx > 0, 



'R a J z|>-_ 

Qafi — o for a ^ /? ± e;, and by the definition of Dj , 

l Rg (x ±Axe t ) 

i=i 



Ax d G /3±e,,/3 = V 7 r- sgn ( 7 f<-) f lg>±A^Q dx ^ q 

,_l V 7 ^-Rsiej AX 



Therefore G^ = G^ + G^' < and G% = G aS + G a ^ > for 

(d) To find the lower bound on G^ we note that J R l Rfi {x+z) — l Rft (x) dx > — Ax d , 
and hence 

The bound then follows since 



/ >Am d v( z )>-f ((M) i w<i(*)+ i i«i>i( j? )) d ^ z )- 



AxG^' >-d[ \z\ dn{z) > -d [ Jd^) 

./4?<UI<i Jo<ui<i -5- 



L 2 

; <|z|<l ^ 0< |^| <1 — 

When J \z\ A 1 d/i(z) < oo, the corresponding bound follows by a similar argument. 



(e) When (jT75j) hold we can estimate Gp t p in the following way 

\z\ c\dz 

Ax \„\d+X 
-<\z\<l — \ z \ J\z\>l 

^(l-m'l+C forA^l, 



G /3,/3>~ /. AlOrfTA - / d M Z ) 



-c A7 ^a d ln4^+G forA=l. 



10 



S. CIFANI AND E. R. JAKOBSEN 



The last equality can be proved using polar coordinates, and o~d is the surface area 
of the unit sphere in R d . Similarly we find that 

- ;*<l»l<l" M"+ A \-a d lnif fo,A = l, 

and since ^1 — (^p) 1- *! is less than 1 or (-^) 1_A when A < 1 or A > 1 respectively 
(and when Ax < 2), the proof is complete. □ 

From the two facts that G£ > when a ^ f3 and sgn(u)A(u) = \A(u)\, we now 
immediately get a Kato type inequality for the discrete nonlocal operator Q. 

Lemma 4.2. (Discrete Kato inequality) If {u ai v a } aeZd are two bounded sequences, 
then 

sgn(u a - v a ) ]T G$(A(up) - A(v p )) <J2 G P \ A M ~ A M\ ■ 

/3eZ d /3<EZ d 

From Lemma \4. II it also follows that the explicit method (13.61) and the implicit 
method (I3.5[) are conservative and monotone, at least when the explicit method 
satisfies the following CFL condition: 

(4.3) cLa — tt \ < 1 where a a is defined in fl4.ll) . 

cr M (Ax) 

Here c is defined in Lemma 1431 and La denotes the Lipschitz constant of A. When 
the Levy measure y, also satisfies (jl.5D , we have a weaker CFL condition 

(4.4) cLa — ti; s < 1 where a\ is defined in (14.21) . 

a x (Ax) 

Proposition 4.3 (Conservative monotone schemes). 

(a) The implicit and explicit methods Q3.5p and Q3.6p are conservative, i.e. for an 
^-solution U, 

E^ = E^- 

a a 

(b) The implicit method is monotone, i.e. ifU and V solve (|3.5p . then 

jjn <yn ^ jjn+1 < yn+1 j Qr ^ > Q 

(c) If (I4.3[) (or (14.4)) and (|1.5[) ) holds, then the explicit method (I3.6[) is monotone. 

Remark 4.4. The CFL condition (|4.3[) implies that < C in general (just as 
for the heat equation), and ^ < G when J" |z| A ldy(z) < oo. Condition (|4.3p is 
sufficient for all equations considered in this paper. In real applications however, 
typically (|1.5[) holds, and the superior CFL condition (|4.4p should be used. 

Proof, (a) Sum Q3.5P or Q3.6P over a, change the order of summation, and use 
Lemma |4~T1 (a): 

E ^ n+1 = E ^ + Ai E ^)( E = E ^ 

(c) Let T a [u] =«„ + At X^ez d the right hand side of (I3.6p . By Lemma 

ED (c), > for a ^ /3 and hence 

d U(3 T a [w] > for /3^a. 

Since ^4 non-decreasing and < 0, we use the lower bound on G^ in Lemma I4TT1 
(c) to find that 

du a T a [u} = l + MG a a A'{u a )>l- cL A - A l 

a^Ax) 
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which is positive by the CFL condition (|4.3p . 

(b) The proof is similar to and easier than the proof of (c). □ 

We then turn to checking the consistency of the method, and to do that we write 
G a>fj = G r aS + G a ,p, T and G Q,/3 = G" Ar + G?' for r > where 

G a,P = TTd / / l Rfi (x + z)-l R p(x) dfl(z)dx, 

' R a J 4?<\z\<r 



Ax 1 

1 

Ax d 



G a ,n,r = — 2 I I ^-R e {x + z) - 1r p {x) d(i(z)dx, 

<R a J\z\>r 



for 7 ; ; 2 = / 3«l|*|<i 

y^<| z |< r 

G r'" = ^ill^ r j R Dtl R M dx. 

If r < we set G^ = = G a ^' r . We also define 

G^ ,r = G^x a -\- G a '^ ,r and G^ r — G Q .^ ;r -+- Gy' , 
and note that Lemmas 14.11 and 14.21 obviously still holds with G^' r or G^ r replacing 

Proposition 4.5. (Cell-entropy inequalities) 

(a) If U is a solution of the implicit method (I3.5p . iften, /or all r > and k € R, 



»7 



fc) < fc) + a* ]T g;^ A(fc)) 



(4.5) 

+ Air/([^+ 1 ,fe) J] G£ r A([/* +1 ). 

(7>J Assume the CFL condition (|4.3[) ('or (14.4)) and (| 1 . 5|) ) holds. If U is a solution 
of the explicit method (|3.6p . i/«en ; /or all r > and k € R, 

<T,(uz,k) + At £ G;' r #((/:),#)) 

(4 6) ^ eZd 

+ At ) /(E£+ 1 ,fc) J] G| )r A(L7«). 

Remark 4.6. In the cell-entropy inequality for the explicit method, the ry'-term 
appears in the "wrong" time. In Section^ we will see that this leads to worse error 
estimates for the explicit method than for the implicit method. 

Remark 4.7 (Convergence to entropy solutions). Proposition 14.51 and a standard 
argument show that any G([0, T]; L 1 1 oc (R d ))-convergent sequence of (interpolated) 
solutions uax of (|3.5p or (|3.6p . will converge to an entropy solution of (|l.ip . Wc 
refer to Theorem 3.9 in [32] and Section 4.2 in [TT] for more details. Convergence 
to the entropy solution also follows from the error estimates of Section [6] 

Proof, (a) By (|3.5p we easily see that for any k £ R, 

K +1 V*<t£V* + At l (fe , +oo) (CC +1 ) t»{A{U n+1 )) a , 

A k > US A k + At l(_ 00 , fe) (C/^+ 1 ) &(A(U n+1 )) a . 
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Subtracting and using rj(u, k) = \u — k\ and T}'(u, k) — sgn (u — fc), we find that 

k) < r,(U:,k) + At V '(U2 + \k) C^(A(U n+1 )) a . 
For any r > 0, we use Lemmas 14.11 (a) and 14.21 with G^' r replacing G^ to see that 

/3GZ d 



= v\U^\k) £ G a r (A(u; +1 )~A(k)) 




< G a ^ v (A(u; +1 ),A(k)). 

The cell entropy inequality now follows from writing Gp = r + G a ^ r and using 
the above inequalities. 

(b) By p.6p and monotonicity (Proposition 14.51 (c)) we obtain the following in- 
equalities: For all r > 0, 

U2 +1 vk<u™vk + AtJ2 G T A (uj} V fc) 

/36Z d 

+Aii (fc , +00) (f/r i ) Y Gp.rMUp), 

pez d 

U2 +1 /\k>U™/\k + AtJ2 G T A (Up A fc) 

/3&Z d 

+ Atl ( _ 00 , k) (u: +1 ) Y G% r A{U%). 

I3ei. d 

Since r](A{U), A(k)) = A(UV fc) — A(U A fc), the cell entropy inequality follows from 
subtracting the two inequalities. □ 

5. A PRIORI ESTIMATES, EXISTENCE, AND UNIQUENESS 

In this section we state and prove several a priori estimates for the solutions 
of the numerical methods (|3.5[) and (13.61) . In what follows, we will use different 
interpolants u of the solutions U™ of the schemes. For the implicit method (I3.5[) 
we take 



(5.1) u(x, t) = U2 +1 for all (x,t) G R a x (t n ,t n+1 
while for the explicit method (|3.6[) . 

(5.2) u(x,t)=U£ for all (x,t) G R a X [t n ,t n+1 ). 
We now prove the following a priori estimates for u: 

(5-3) ||M(-,t)|| L i( R d) < ||u || L i( R d), 

(5.4) \\u(-,t)\\ L oo {Rd) < \\u \\ L ^ (Rd) , 

(5-5) \u(-,t)\ BV(Rd) < \u \ B v(WL d )- 



Lemma 5.1. (A priori estimates) 

(a) If U solve (|3.5p and u is defined by (15. ip . then the a priori estimates (|5.3p - 
(153)1 hold for allt>0. 

(b) Assume the CFL condition (g3J (or and ([T3j) ) ZioWs. J/ [/ sofoe jXB} 
and u is defined by (|5.2p . i/ien i/ie a priori estimates (|5.3p - (|5.5p ZioZd /or allt > 0. 
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Proof. Since the schemes are conservative and monotone, cf. Proposition 14.31 this 
is a standard result that essentially follows from the Crandall- Tartar Lemma. For 
explicit methods in part (b) we refer to e.g. Theorem 3.6 in [22] for the details. 

We did not find a reference for implicit methods, so we give a proof of part (a) 
here. See also [T7j for the case when A is linear. Let u a = U™ +1 , h a — U™, and 
write (13.51) as 

(5.6) u a - At J2 G ?3 A (up) = h a . 

/3ez d 

We prove (|5.3p . Multiply (|5.6[) by sgn(uo,) and use Lemma l4~2l to get 

|u„|- At J2 G?\Mup)\ < K\> 

i3ei d 

which by Fubini's theorem and the fact that J2aei, d G % ~ ^ implies that 

53 E i^i- 

By the definition of u a , h a and an iteration in n, it follows that 

E ™ ^ E i^i- 

aSZ d a£Z d 

By (EH), \\u(;t)\\ Ll{Rd) = Ax d J2 a&d \US\ for t € (t n ,t n+1 ], and (53} follows. 
To prove (|5.5|) . we subtract two equations (|5.6[) evaluated at different points, 

u a - u a - ei - At Y ( G A ( u p) ~ G a f ei A(upf) =h a - h a - e , 

i3ei d 

and use the fact that G2 = Ga+tl *° see 

U a - U a - ei - At J3 G /3 ( A ( u p) - M u P-eif) =h a - h a - er 

fiez d 

Then we multiply by sgn(u a — u a - ei ), use Lemma |4"T21 and sum over a, to find that 

E l U " - U Q _ ei | < >^ \h a - h a - ei I ■ 
aGZ d a<EZ d 

The estimate (|5.5p then follows by iteration and the definitions of u a ,h a , u. 

It remains to prove (I5.4p . Note that since ^ a \ u a\ < °o by (|5.3p . there is an 
«o such that sup Q u a — u ao . Moreover, the parabolic term is nonpositive at the 
maximum point: since X^ez d G % = ® an d S,3ez d l^sl < 00 > 

£ A( U/j ) = ^ G?(A{up) A{u ao )) < 0. 

,8£Z d i3ez d 

Then by the above inequality and (|5.6I) . 

sup u Q = w Qo < Uqq - At 53 G< j3° A ( u fi) = h a < sup h a . 

aSZ d p £Zd a£Z d 

In a similar way we find that inf ag zd h a < inf ag g<i w a and (I5.4p follow from the 
definitions of u ai h ai u and an iteration in n. □ 

Lemma 5.2 (Global existence and uniqueness). 

(a) There exists a unique solution U" € I 1 of the implicit scheme (|3.5I) for all n > 0. 

f&j Assume the CFL condition (14.31) (or (|4.4p and (|1.5p ) holds. Then there exists 
a unique solution U n £ l l of the explicit scheme (|3.6p /or all n > 0. 
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Note that U n £ I 1 implies that u(-,t) £ L 1 (R d ). 
Proof, (a) Let u a = and h a = U2, rewrite {33} as (HHJ), define 

T Q [u] = u Q - e( u a - At ^2 G< p M u p) ~ h a 

^ /3ez d 

and let e be such that 

e(l+L A c , -j ) < 1. 
V cr M (Aa;)y 

We first show that T Q is monotone, i.e. u < v implies T a [u] < T a [v]. For a^(?, 
Gg > by Lemma B~T1 and hence since A non-decreasing, 

d U pT a [u] >0. 

Moreover, since A non-decreasing and — g— £ < 0, 

At 



d Ua T a [u] = l-e + eAtGZA'(u a )>l-e[l+L A c^ 

which is positive by our choice of e. 

Since T is monotone and A is nondecreasing, 

(r a [u] - T a [v\Y < v«] - r a [«]) 

a a 

= (1 - e) ^KV!) a -jJ a ) + £A^ G ? ( A ( w V V P) ~ A{ - v ^ 

a£Z d aeZ d /3GZ d 

= (i - e) - + e A< E ( E G i) ( a m - ^m) + - 

a£Z d /3eZ d aGZ d 

A similar estimate holds for $Z a (T Q [u] — T a [u])~, and since J2 a ez d G p ~ ^> we have 
shown that 

E l T «N " T «MI < (1 - e) E l M « " 

ctSZ d a£Z d 

So T Q is an ^-contraction and Banach's fixed point theorem then implies that there 
exists a unique solution u £ I of T a [u) = u a and hence also of (15. 6p . 



(b) Existence follows by construction and the a priori estimates in Lemma 15.11 
Uniqueness essentially follows by monotonicity and J2 a g b — ; Assume two so- 
lutions U n and V n , subtract the two equations and multiply by sgn([/ n — V n ), 
and use the Kato inequality (Lemma 14. 2[) along with J2 a G< p = to show that 

E a P n -v n \ <E a \u°~v°\. ^ □ 

We have the following regularity estimate in time: 
Lemma 5.3. (Regularity in time) 

(a) Assume (A. 2) - (A. 4) hold, and let U be a solution of the implicit method (|3.5[) 
and u defined by (|5.ip . Then 

\\u(; S) - u(-,t)|| il(R d) < CT M (|S - t\ + At) 

for all s, t > 0, where 

I */Jm >0 M A1 d M(>) < °°; 
<W) = \ 

fr otherwise. 
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(b) Assume (A. 2) - (A. 4) and (|4.3|l (or and (|T3)) ) hold, and let U be a 

solution of the explicit method (|3.6[) and u defined by (15. 2[) . TTien 

||u(-,s) - u(-,t)\\ L i {Rd) < ct p (|s - t\ + At) 

for all s,t>0, where is defined in (a). 

Proof. The two proofs are essentially identical, so we only do the proof for case (a). 
1) By (|3.5p . we find that for any x € R a , 

At r 

u: ~ c/r 1 = j R c[A(u n )]( X )dx. 

Take a test function < (f> <E and define <^> Q = -^p J fl <Ky)dy ar± d </>(a0 = 
^ a </)Q,lfl a (x). Multiply the equation by Ax d (f> a and sum over a to find that 



<P(x)(U n (x) - ET^fc)) dx = At / 0(a;)£[A(C/™)](x) dx, 
where Let C* be the adjoint of C, then since U is constant over R a , 

(x)(U n (x) - U n -\x)) dx= f 4>(x)(U n (x) - U n -\x)) dx 
= Atf (4>{x)-<j>(x))C[A{U n )]{x)dx + At( C*[<jy}(x)A(U n )(x)dx. 



2) Let u> e be an approximate unit, i.e. uj e (x) = j^ w (f ) where < oj € and 
j &d uj dx = 1. Take = u; e (y — x) in the equation above and let U™ — U n * cu £ : 

ftn _ fjn-i = A ^- £ _ ^ £[A(U-")] + Ai£*[u; £ ] * A(U n ). 

By Fubini we then find that 

^||L> - ^lUi <\\Q e - u e \\ L 4£[A(U n )]\\ L i + \\£*[ Ue ] * A(U n )\\ L1 =h + I 2 . 

3) To estimate Ji, note that by a standard argument 

||u> £ - u E \\ L i < \uj £ \bvAx = —Ax, 

£ 

and then by the definition of £ in f|3 . 3[) . Fubini, the L 1 n BV regularity of U n 
(Xemma l5.il) . and the regularity of A in (A. 2), 

\\C[A(U n )]\\ L1 

[ A(U n (x + z)) - A(U n (x)) - z ■ D Ax A(U n (x))l lzl<1 dxdfi(z) 

>^f- JR d 

< f (2\A(U n )\ BV \z\l lzl< i+2\\A(U n )\\ L il lz>1 )dti(z) 
<C I b|Ald/i(z)<— / \z\ 2 A 1 d/i(z). 

J\z\>£f Ax J\z\>0 

These estimates along with (A. 3) shows that I\ < Ce^ 1 . 
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4) Then we estimate I2. Note first that since Dax = D + (Dax — D), we can use 
Taylor's formula to see that 

4>{x + z) — <f>(x) - zD& x 4>{x) 

= I (1 — s)z T D 2 cf>(x + sz)z ds ± Ax Zj I (1 — s)(p XiXi (x ± sAx) ds. 
Jo i=1 Jo 

This identity along with the definition of C* , repeated use of Fubini, and one 
integration by parts in x, then leads to 

t*[Lo e ] *A(U n )(x) 

(1 - s)Duo e {x - y + sz)z ® z DA(U n (y)) dydfi(z) ds 



J^<z<l 



d r i r 



'<|z|<l 



z >1 



(1 - s)d Xi w s (x -y± sAx)zi d Xt A(U n (y)) dydfi(z)ds 
(u e (x -y + z)-u> e (x- y))A(U n (x)) dfi(z) dy. 



Here DA(U n (y))dy should be interpreted as a measure, and J \DA(U n (y))\ dy = 
J d\A(U n )\(y) = \A(U n )\BV- By Young's inequality for convolutions (Fubini in our 
case), we then find that 

h < %e\Bv\A{U n )\ BV I |z| 2 d^) + 2||^|| ii p({7")|| L1 f dfl(z) 

z J0<\z\<l J\z\>l 

< Ce-\ 

Here again we have used the properties and regularity of /z, A, U n , and lj £ . 

5) By steps 2) - 4) we can conclude that 

n ^ 
\\Ul l -Ur\\^< £ m-Ut l \W<-\n-m\At, 

j= m+l 

where the constant C does not depend on n or m. By the triangle inequality and 
standard BF-estimates, it then follows that 

\\u n - U m \\ L i <\\U n - U?\\ L i + \\U? - Ul n \\ L i + \\u™ - U m \\ L i 
< \U n \ B ve + -|n - m\At + \U m \ BV e, 

£ 



and hence by taking e = C \/\n — m\ At, 



\\U n ~ [7 m || L i < CyJ\n-m\bt. 

For the time- interpolated function u defined in (|5.1[) . wc then find the following 
estimate 



\\u(;t)-u(;s)\\ L i = \\U n -U m \\ L i < CVjn - m\At < Cy/\t -s\ + At. 

The equality follows since for each t, s there are n, m such that u(x, t) = U n {x) and 
u(x, s) = U m (x). Moreover, by the definition of u, \n — m\ At < \t — s\ + At. 

It remains to prove a better estimate for the case when J \z\ A 1 d/i(z) < oo. This 
proof is similar but much easier than the proof above, so we skip it. □ 
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The time regularity result in Lemma 15.31 is not optimal for Levy operators C 
with order in the interval [1,2). To get optimal results we need more detailed 
information on the Levy measure /i than merely assumption (A. 3). We will now 
prove an improved time regularity result for fractional measures (jl.5p . In this result 
we will need the following CFL condition, 

(5.7) C A^7X< 1 fOT Ae (°' 2 )- 

Lemma 5.4. (Time regularity for fractional measures) Assume (A. 2) - (A. 4), and 
HE5H hold. 

(a) If the CFL condition (|5.7[) hold, U is a solution of the implicit method (|3.5[) 
and u its interpolation defined by (15. ip . then for all s,t> 0, 

\\u(-,s)-u(;t)\\ L1{WLd) <a x (\s-t\+At); <j x {t) = 

(b ) If the CFL condition (14.41) hold, U is a solution of the explicit method (|3.6p and 
u its interpolation defined by (15.21) . then for all s,t > 0, 

( t A<1, 
||w(-,s)-u(-,t)|| L i( E<J ) < ax(\s-t\+At); cr A (r) = < r a for any a e (0, 1) A = 1, 

[ t* A > 1. 

Note well that in this result we need the CLF condition also for the implicit 
scheme. The reason is that the time-regularity is linked through the equation to 
the approximate Ax-depending diffusion term as will be seen from the proof. For 
the implicit scheme, we can have better results for A = 1 since we can use the less 
restrictive CFL condition (|5.7[) . 




Proof. The result for A < 1 is a corollary to Lemma [5.31 The proof for A > 1 is 
the same as the proof of Lemma 15.31 except that we use different estimates for I\ 
and Ii in step 2). From step 3) in that proof and (| 1 . 5|) and a simple computation 
in polar coordinates, we get that 

^Ax f -i i / \ ( f \z\dz „ 

h<C—J \z\M&n(z)<C— (J a tt r 



e 

~ |~|^ 2 

„2-A 



^ C I As + Ax 2 - A for A > 1, 
_ £ 1 Aa; - Ax In Ax for A = 1. 

To estimate I2, we use Taylor expansions and integration by parts to find that 

C*[u) e ]*A(U n )(x) = 



^^<|z|<e 
1 



(1 - s)Duj s (x -y + sz)z ® zDA(U n {y))dydn(z) ds 

1 

I (u e (x -y + sz)- uj e (x ~ y))zDA{U n (y)) dyd^(z) ds 
:i Js. d v ' 

TAxJ2 / (l-s)d Xl uj e (x-y±sAx)z l d Xt A(U n (y))dydfi(z)ds 
tT] f f [ (u £ (x-y±sAx)-Lj e (x-y))z t d Xl A(U n (y))dydfi(z)ds 

J—1 Jo Je<\z\<l JR d V ' 

(u £ (x -y + z )- lo £ {x - y)^jA(U n {x)) d/J,(z) dy. 



J e< z <1 



z >1 
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Then by Fubini, the definition of u £ , and the change of variables (x, z) — > (ex, ez), 

\z\ 2 dz 



\Dlj s (x + sz)\\z\ 2 dfi(z) < c A e 1 - A / \Du\ dx . 

E <|z|<£ JR d J0<\z\<l \ z \ 

By similar estimates and Young's inequality for convolutions we find that 

+ 2||i4(^)|Ui||w e |Ui / d/*(«) 
J| z |>i 

f £ r 1 - A -1-1 X 1 

< C 





Note that the lne-term comes from the integral over 1 < \z\ < ~. 
As in step 5) in the proof of Lemma 15.31 we then find that 

\\U n - U m \\ L i < \U n \ BV e +\n- m\At(h + h) + \U m \ BV e. 

To conclude, we assume that Ax < e which means in particular that 

h + h < C 

When A > 1, the final result follows from taking e — c(\n — m\At)i and arguing 
as in the end of the proof of Lemma 15.31 Note that in view of the CFL conditions 
(|4.4p and (|5.7|l . the constant c can be chosen such that Ax < e. For A = 1, we 
can use e = c\n — m\At for the implicit method in view of (|5.7[) . and by (|4.4[) . 
e — c(|n — t?t,| Ai)" for any a' £ (0, 1), will do the job for the explicit method. □ 

By the a priori estimates Lemma 15.11 and 15.31 and Kolmogorov's compactness 
theorem (cf. e.g. [221 Theorem 3.8]), we find subsequences of both methods (|3.5p 
and (|3.6[) converging to some function u. The function u inherits all the a priori 
estimates of u, and it will be the unique entropy solution of (jl.ip by Remark 14.71 
In short, we have the following result: 

Theorem 5.5. (Compactness) Assume (A. 2) - (A. 4) hold. If either 

(i) U is the solution of the implicit method (|3.5[) and u defined by (|5.1[) . or 

(ii) U is the solution of the explicit method (|3.6[) , u defined by (|5.2I) , and (|4.3[) 
(or (l4~4f and l|1.5|) ) also holds, 

then there is a subsequence of {u}/\x>o converging in C([0, T]; L 1 (R d )) to the unique 
entropy solution u of (|1.1|> as Ax — > 0. Moreover, 

u e L°°{Q T ) n C([0, T];L 1 (R d )) n L°°(0, T; By(M d )). 

Remark 5.6. This result provides a proof for the existence result Theorem 5.3 in 
[TT] for L 1 l~l L°° n entropy solutions of (jl.ip . and then the general existence 
result in L 1 n L 00 follows by a density argument using the ^-contraction. 



6. Error estimates 



In this section we give different error estimates and convergence results for our 
schemes, estimates that are valid for general Levy measures and better estimates 
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that holds for fractional measures satisfying (|1.5[) . To give the general result, we 
need the following quantities: 

ir = - 1 \z\ 2 Mz), 

6 J\z\<r 

'*=(V + f)(/ <ws J« |d "« + / s|>1 <""4 

I3 = £At(u) / d/x(z). 

J\z\>r 

Theorem 6.1. (Error estimates) Assume (A. 2) - (A. 4) hold, and let u be the 
entropy solution of (jl.ip . 

(a) Let U be a solution of the implicit method f|3.5[) and u defined by (|5.ip . Then 
for all e > 0, < 5 < T, and ^ < r < 1, 

(6.1) \\u(-,T) - u(; T)\\ Ll{Rd) <C T (e + £ g (u) V £ 5 {u) + I{' r + I'/' r 

(b) Assume also (|4.3I) holds, and let U be a solution of the explicit method (|3.6[) 
and u defined by (|5.2[) . Then for all e > 0, < S < T, and < r < 1, 

(6.2) |K-,T) - T)|| il(Rd) < C T ^ + f ff (ti) V S s (u) + I{' r + I £ 2 5 ' r + I^j , 

The proof of this result will be given in Section [5] 

Corollary 6.2 (Convergence). Under the assumptions of Theorem \6.1\ the solu- 
tions of the implicit method (|3.5p and the explicit method (|3.6[) both converge to the 
unique entropy solution of p. ID as Ax, At — > 0. 

Proof. The result follows from the error estimates of Theorem 16.11 by first sending 
Ax, At — > 0, then r — > 0, and finally e, S — > 0. □ 

We will now see how Theorem 16.11 (along with Lemma 15. 4[) can be used to pro- 
duce explicit rates of convergence for our scheme in the case of fractional measures 
satisfying (jl.5p . First we define 

' r\ Ae (0,1), 

(6.3) o{ U {r)={ r^|logr| A = 1, 



2-A 



and 



(6.4) a f x {r) = { T La 

\ T** Ae (|,1)U(1,2). 

Theorem 6.3. (Convergence rate for fractional measures) Under the assumptions 
of Lemma \5.4\ ( including (jl.5p and a CFL condition for the implicit scheme ), for 
all A G (0,2), 

{Ct cr{ M (Ax) for the implicit method (|3.5p , 
Ct fx ( f or the explicit method (|3.6D . 

Note that the rate for the explicit method is worse due to the extra term ig in 
Theorem 16.11 



t— AG (1,2), 
AG(0,|], 
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Corollary 6.4 (Explicit scheme when A = 1). Let the assumptions of Lemma \5.4\ 
(b) hold with A = 1 and let a 6 (1,2) be arbitrary. If the stronger CLF condition 
C-^ < 1 holds, then 

\\u(-, T) — u(-, T)\\ L imd) < Ct cr® x (Ax) for the explicit method (|3.6D . 



Proof. Note that the CFL condition (|4.4[) is satisfied and that the assumption (jl.5p 
holds with any A 6 [1,2). Hence the result follows from the A > 1 case in Theorem 

Efl □ 



Proof of Theorem 1 6'. ffl Let us first give the proof for the implicit method (|3.5 
First we note that by Q1.5II . 



\z\ 2 dfi(z) <c x [ dz<0 (r 2 - x ) for all A e (0, 2), r < 1, 

>\z\<r J\z\<r \ z \ 

while 

'0(1) if A 6 (0,1) 

\z\ dfj,(z) < c x / TZHTx dz = < 



i-<|z|<1 Jr<\z\<l 



0(\lnr\) if A = 1, 
0(r 1 - A ) if A e (1,2) 



Using these estimates along with the CFL condition (|5.7[) and Lemma 15.41 we find 
that the estimate (|6.1[) in Theorem 16.11 takes the form 

- e + (5+ rf_i + (^ + ^) if A €(0,1), 

IK, T) - u(;T)\\ LHm < C T { e + 8 | hxS\ + f + | lnr| + ^) if A = 1, 

e + §i + t^ +r l-X^ + A^ if A 6(1,2). 

The conclusion then follows by taking r = Ax for all A 6 (0, 2), e = 8 = \/ Ax for 
A 6 (0, 1], while e = Ax 2 ^ and 8 = Ax? for A 6 (1, 2). 

For the explicit method (|3.6|) we also need to take into account the extra /3-term, 



I%=o- x (At)[ dfi(z) 

J \z\>r 



>r 
0(r- 



Lemma r5.4[ and the slightly more restrictive CFL condition (|4.4[) . The expression 
(|6.2[) in Theorem 16.11 then takes the form 



\\U(;T)-U(;T)\\ L1 

\ _L ( Ax _j_ _ 

S 

f: + 5i + rL± +r i-A^ + ^ + ^ if A 6 (1, 2 



< C T 



e + 5+ r — + (^+£f) + Az if A e (0,1), 



We minimize two and two terms and take the maximum minimizers, first w.r.t. e 
and 8 and then w.r.t. r, 



\\u(-,T) - u(-,T)\\ L i (VLd) 

V 2 ^ +Ax? + if A 6 (0,1) 

r ¥ + r^Ai? + r T+iAiW + ^ if A e (1,2), 
Aa;5 + Az^tx if A 6 (0, 1), 

Aa; 2 ^ + As" + As^ if A 6 (1,2). 



< C T 



< C T 
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The final result follows since §5^ > §=^ for A 6 (^,2) and < \ for A e 
(|,2). " 1 ~ ~ " " □ 

Remark 6.5. The rates can not be improved by taking a different truncation of the 
singularity, i.e. replacing in the method 

G ° s >l 3 = ~^dJ J A ^-B^{x + z) - l Rt) {x) A^{z)dx 

by 

G at p = -r—7 / 1^ (a; + 2) - l Rfl (x) dfi(z) dx. 

The reason is that the function p\ that minimize the error expression 



Pa" A (Ax) 1-x ,/AiA 



+ Pa - A (^) 



l.r 



<5 



is always p\(Ax) = 0(Ax)\ 

Remark 6.6. We believe that the rates for the implicit schemes are optimal, at 
least when there are nonlinear convection terms in the equation (i.e. when / ^ in 
(II. ip . see Section [7]). But we have not found analytical examples confirming this, 
nor have we been able to observe the above rates in preliminary, but probably too 
crude, numerical tests. Maybe it is not straight forward to construct analytical 
or numerical examples confirming the optimality of the rates. We leave it as a 
challenge for people with more experience in realizing numerical schemes to test 
the optimality numerically. 

7. Convection-diffusion equations 

In this section we discuss how to extend the results established in the previous 
sections to the case / ^ 0. Note that all the arguments needed to handle the 
additional /-term are well-known. We consider the following numerical methods 

d 

(7.1) U» +1 =K + At ^ D T fi > u atli ) + At^(A(U n+1 )) a , (implicit) 

i=i 

d 

(7.2) U2 +1 = K + At .Df fi(U£, f/£+e, ) + At&(A(U n+l )) a , (expl-impl) 

i=i 
d 

(7.3) U2 +1 =K + At^DTfiiK, U: +ei ) + Atfr(A(U n )) a , (explicit) 

i=i 

where 

(i) D^U a = -£-(U a — U a - ei ) and {e{\i is the standard basis of R d , and 

(ii) f — (fi, . . . , fd) is a consistent (i.e. f(u,u) = f{u)), Lipschitz continuous 
numerical flux which is non-decreasing w.r.t. the first variable and non- 
increasing w.r.t. the second one. 

Remark 7.1. Some examples of numerical fluxes / satisfying (ii) are the well-known 
Lax-Friedrichs flux, the Godunov flux, and the Engquist-Osher flux, cf. e.g. [22]. 

For the schemes (|7.2D and (17.31) . we also need the CFL conditions 

(7.4) 2dL F ^- + cL A - - s < 1 and 2dL F —^<l 

Ax a^(Ax) Ax 

respectively (compare with (|4.3p ). where is defined in (|4.1[) and Lp is the Lips- 
chitz constant of /. Then the all the a priori estimates and other results of Section 
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[5] continue to hold for the new schemes, and we still have compactness via Kol- 
mogorov's theorem. The modifications needed to identify the any limit as the 
unique entropy solution of are standard and can be found e.g. in Chapter 3 
in [22], and hence the convergence of the methods (|7.ip - (|7.3p follows. 

We will now give the statement of the result of Theorem 16.11 that is valid for the 
current setting where / 7^ 0. To do so we reuse the quantities I e { r and /J of section 
[6j but redefine I^ &,T as follows 



= I — + -T" ) ( 1 + / \Z\ d//(z) + / d/i(z) 

V 6 J \ Jr<\z\<\ ->\A>^ 

Theorem 7.2. (Error estimates) Assume (A.l) - (A. 4) hold, and let u be the 
entropy solution of 

(a) Let U be a solution of (|7.ip or (|7.2[) and u defined by (|5.ip . For (|7.2p we 
also need the second CLF condition in (|7.4p . Then for alle>0,0<5<T, and 
*f<r<l, 

\\u(-,T) - u(-,T)\\ LHRd) <C T [e + S 5 (u) V £ g (u) + I{' r + I e /' r ^j . 

(b) Assume also that the first CFL condition in (|7.4p holds, and let U be a solution 
of (EH]) and u defined by (g^). Then for all e > 0, < S < T , and < r < 1, 

|K, T) - u(-,T)\\ LHRd) <C T [e + £ s (u) V £ s (u) + F{ r + I £ /' r + L 

The proof is essentially equal to the proof of Theorem l6 . 1 l augmented by standard 
Kuznetsov type computations to handle the /-term, cf. e.g. |HJ Example 3.14]. 
We skip it. 

Remark 7.3. It is easy to see that the contribution to the error from the discretiza- 
tion of the /-term is always less or of the same order as the contributions of the 
other terms. In particular, for fractional measures (jl.5p . we immediately get that 
the schemes satisfy the error estimate of Theorem 16.31 with modulus a{ M for (|7.ip 
and (j7T2j) and modulus of x for (|7T3|) . 

8. The proof of Theorem 16.11 
Proof of Theorem \6.1\ for the implicit method (|3.5p . 

1. We use Lemma l2~2"l to compare the solution of the scheme to the exact solution. 
In the resulting inequality, we introduce the scheme via the time derivative and the 
initial/final terms. To do this, we use integration by parts on each interval {t n , t n+ i) 
and summation by parts to get discrete time derivatives on u so that we can use 
the cell entropy inequality (|4.5I) . We get that (remember the definition of u) 

•n(u(x, t), u{y, s)) dt<p e ' S (x, y, t, s)dw + initial and final terms 

JV-l 



Aa; At 



f EE Uu- + \u(y,s))- V (U-,u(y,s))) [ 

JJqt n=0 Q6Zd v / JRo 



ip e ' S (x, y, t n+ i , s) dx dyds. 



Let (p e ' S = ip e ' S (x, y, t, s) be the function which for each (y, s) € Qt is defined by 



— / (p e ' 5 (x,y,t n ,s) dx for x G R a , t G (t n -i,t n ], 
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and use above equation along with the cell entropy inequality (|4.5p and Lemma |3. II 
to write the inequality of Lemma 12.21 in the following way 

\\u(;T)-u(-,T)\\ L i m < C T (Ax + e + £ s {u) V £«(«)) 

r)(A(u(x, i)), A(u(y, a))) Cf [^ s {x, -,t, s)](y) dw 



Qt jjq 



Hi 

r)(A(u(x, t)), A(«(y, «))) Cf [<^(-, y, t, s)](z) dw 

Qt JJQt 

J 

V 

ff 2 

T/(«(a:, t), u(y, s)) £"' r [A(u(-, t))](z) (^' 5 - ^' 5 )(x, y, t, s) dw 

Qt JJQt 



H 3 

,<5 V7 ,^,5\ 



v(A(u(x, t)),A(u(y, a))) ^ > r ■ (D^ 6 - V x ^ A ){x, y, t, s) dw . 

Qt 



Here we have also used the notation 

C[<j>](x) = CM(x) + + 7^ • D Ax cf>(x) 

where C is defined in (|3.3p . C r — C for r > 4r, and 



= / ^(ac + 2r) - 0(a;) - l| z i<iz • D& x <t>(x) dfi(z). 

J^<\z\<r 

Note that the discrete operator £>i = Dax,i (see p. 21) ) always acts on the x-variable 
(the variable of u). To complete the proof we need to estimate Hi, . . . ,H±. 

2. Estimates of Hi and Hi- By Taylor's formula with integral remainder, inte- 
gration by parts, and Fubini (- see e.g. Lemma B.l in [2j for more details), 

\Hi\<[( I! I I (l-T)\D yV (A(u(x,t)),A(u(y,s)))\ 



T J\z\<r JO 



us(t — s) \D y £l e (x — y + tz) | \z\ 2 dr d(j,(z) di 



< \h A I \u(-,s)\ BV ma)ds I \D x Q e (x)\ dx I u E (t)dt I \z\ 2 d(i(z)) 



z\<r 



|_D x n e (x-y+rz)| 

^La f \v,(-.s)\ n „ m ^ds [ 
< C T L A \u \ BV{Md) / |z| 2 d^(z). 

J\z\<r 

Here we also used Theorem l2.1l and the standard estimate J Rd \D x Q e (x)\ dx = O(-j). 

We find a similar estimate for Hi via a regularization procedure and the ar- 
gument for Hi above. Let (p g ' S be a mollification in the x- variable of (p e ' S , i.e. 
(p e g ' S — ((f^ * x Q p where the convolution is in x only. Then (p g ' S is smooth in x, and 

\^ e g S (-,y,t,s)\ BV(Rd) < \^' S (-,y 7 t,s)\ BV{Rd) < \ip e > s {- 7 y,t 7 s)\ BV{Rd) = O (e _1 ) , 

where the first inequality holds for all g small enough (cf. e.g. [301 Theorem 5.3.1]), 
while the second one is obvious. Let us call 

H$= H II v(A{u{x,t)),A(u(y,s)))Cf[^/(;y,t,s)}{x)dw. 

JJ Qt J J Qt 



24 



S. CIFANI AND E. R. JAKOBSEN 



First note that lim e _ i .o i?| — H2 by the dominated convergence theorem since we are 
integrating away from the singularity and (p e g ' S (• ,y,t, s) — > (p e ' S (• ,y,t, s) pointwise. 
Then, since (p 6 ' s (-, y, t, s) is smooth, we repeat the argument used for Hi and obtain 



\H§\ < C T L A \uo\ BV(] g^ \ip e < \ BV md) / \z\ dfi(z). 

J*f<\z\<r 

Since \tp e g ' S \BV(K d ) = ^(e -1 ), we can take the limit g — > and get 
\H 2 \ < C T L a |u |w(R d ) / kl 2 dfi(z). 

J\z\<r 

3. Estimate of H3 . By the definition of (p e ' S and properties of mollifiers, a standard 
argument shows that 



ip e ' (x,y,t,s) - ip e ' (x,y,t,s)\ dyds 



/Ax Ai\ 

< d\n e \ BV \\ujs\\ L iAx + d\\n e \\ L i\cjs\BvAt < o f — + — J . 

Similar estimates are given in e.g. |13j . This estimate along with several applica- 
tions of Fubini's theorem then show that for all < r < 1, 

|if 3 |< ff \C^ r [A(u(;t))](x)\ t ff \^ s (x,y,t,s)-^ 5 (x,y,t,s)\ dyds ) dxdt 



< cLa ( 1 — t- ) [ II I \u(x + z,t) — u(x, t)\ dfj,Cz) dx dt 

• e 5 J \ JJ Qt J r< \ z \<i 



\ 



( \u(x + z, t) — u(x, t)\ d/j,(z) dx dt 

Iq t J\z\>1 J 



< C T L A (— + -r) f \ u o\bv f \z\dfi(z) + \\u \\ L i f d[i(z 

V 6 \ Jr<\z\<l J\x\>l 



4. Estimate of H4. Let I £ (0, . . . , d) and write 

N- 1 rtn 



Hi,i = lf' r ff EE v(A(UZ),A(u(y,s))) f b^&y^s) dw 

JJOt n Jt„ JR„ 



H 

N-l 



tf V ff EE r,(A(U2),A(u(y,s))) f " +± f d Xl ^ s (x,y,t,s) dw . 

JJ Qt aeZd n =0 ■ / *« jR « 



H 2 

A 1 



Since J t n+1 f R (p € ' S (x,y,t,s)dxdt = J t n+1 J R <// :<5 (x, y, tn+i, s) dxdt by definition, 
we can use summation by parts to find that 



N— 1 t 

H h =~ II EE A*7(A(I£), *))) f " +1 f ^(x, y, t n+1: s) dw. 
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2o 



Integration in the a^-direction followed by summation by parts leads to 



N-l 



Hi = -Ax //EE Div(A(US), A(u(y, a))) 

JJQt adZ d n=0 



(p e ' S (x\ Xl = x , y, t, s) dx%... dx/_i dx ;+ i . . . dx d dt dy ds. 



Here we first integrated d Xl (p e ' S (-, y, t, s) along the interval (x ai , x a[+1 ) to obtain 
the difference (p e ' 5 (x\ xl — Xai j , y, t, s) — ip € ' S (x\ Xl — Xai , y, t, s), and then we used sum- 
mation by parts to move this difference onto r](A(U%), A(u(y, s))). Note that 
x \xi=x a = ( x i-> ■ • ■ j x l-l> x a t , xi+i, ■ ■ ■ , Xd), and that xi — x ai is fixed here while 
the other variables Xj ,3^1 vary. 

By the above computations, the inequality \Dir}(A(U%), A{u{y, s)))| < \L>iA(U™)\ 
(i.e. | jet — k\ — \b — k\\ < \a — b\), and Fubini, we find that 



N-l 



V ff E £ A»?(>l(E£M(«(lM))) 

JJQt a£Z d n=0 

(tp e,s (x\ xl=Xai ,y,t,s) - tp e > s (x a ,y,t n+1 ,s)\ dxdt dyds 



N-l 



< if'* E E \DiA(K)\ 



W ,5 { x \xi=x a , , V, t, s) - ip e ' s {x a ,y, t n+1 ,s) 



R a JJQ 



dw. 



Since (f> e ' S (x, y, t, s) — Q e (x — y)ojg{t — s) and (x, t) G R a x (t n , t n +i], we find as in 
part 3 that 



Qt 



^(xix^xc^y,*^) - ip e ' d (x a ,y,t n+1 ,s) 



C,S/ 



dyds 



\ ( Ax 


At\ 


) =0 (-" 


--) 







Summing over I we then find that 
\H4<dC\ 7 »*> r \ 



^ + ^)(EEiA^raiAtA^), 

' V n=0 a£Z d ' 



and since J2 a ez d \DiA(U™)\Ax d = \A{u{- ,t n ))\ B v < L A \u \ B v, we conclude that 

\z\dfi(z). 

In view of part I - 4 the proof is now complete. □ 



|/./,| C,L A ( ^ + ^ 



r<|*|<l 



Proof of Theorem \6. 1\ for the explicit method (|3.6[) . We argue as in the beginning 
of the proof for the implicit method, replacing the implicit cell entropy inequality 
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by the explicit one (|4.6[) , and find that 

\\u(-,T) - u(-,T)\\ L1(Rd) <C T (Ax + e + £ s {u) V £ s (v)) 

v{A(u(x, t)),A(u(y, a))) [<p*\x, ■, t, s)]{y) dw 



Qt JJQt 



+ // r,(A(u(x,t)),A(u(y,s)))Cf [^ d (;y,t,s)](x)dw 

JJQt JJ Q t 

+ !! [[ r/(u(x,t + At),u(y,s))C^ r [A(u(;t))](x)^ 5 (x,y,t, S )dw 

JJ Q T JJ Q t 

- [I [I r ] '(u(x,t),u(y, S ))C^[A(u(;t))](x)^ s (x,y,t, S )dw 

JJ Qt JJ Qt 

V (A(u(x, t)),A(u(y, a))) 7 A '*' r • - V x ^' & ){x, y, t, s) dw. 



I Qt JJ Qt 

The difference with the previous proof is the interpolation (I5.2j) , and more impor- 
tantly, the new £ M ' r -terms. Note that by a change of variables, 

ff [[ V '(u(x,t + At),u(y, S ))£^ r [A(u(;t))](x)^ s (x,y,t, S )dw 

JJ Q T JJ Qt 

= ft [[ v'(n(x,t), U (y,s))/Z^ r [A(u(;t-At))}(x) ( p e ' 5 (x,y,t,s)dw 

JJ Qt JJ Qa±,t 
T+At 

r)'(u(x, t),u(y, s))£^ r [A(u(-,t - At))] (as) & e ' 5 (x, y, t, s) dw, 
where Q a ,b — K d x (a, b). The last term on the right can be estimated by 
At\\4> s ' 5 \\ L i l\A(u)\ BV [ \z\dn(z) + 2\\A(u)\\ L i f dfi(z)) 

\ Jr<\z\<l J\z\>l J 

= 0{At) [ \z\dn(z). 

Jr<\z\<l 

By similar computations, we can write the £ M ' r -terms in the above inequality as 

ri'(u(x, t), u(y, s)) C^ r [A(u(-,t - At)) - A(u(; t))](x) <?> s (x, y, t, s) dw 



/ 



V '(u(x, t),u(y, s)) C^ r [A(u(- 7 t))](x) - ^ s )(x, y, t, s) dw 

IQt JJQt 
+ 0(At) [ \z\dfi(z). 

Jr<\z\<l 

Here we estimate the first term using the time regularity of u, 

I<[[ \C»> r [A(u(;t-At))-A(u(;t))](x)\ [[ ^'\x,y,t, s )dw 

JJ Qt JJ Qt 

V * ' 

=0(1) 

<C2L A ( [ \\u(;t-At)-u(;t)\\ L i m dt][ dp(z) 
Wo / J\z\>r 



<C T £ M {u) / dfx(z), 

J\z\>r 

where £At(u) is defined in (I2.2j) . Now all the remaining terms can be estimated as 
in the proof for the implicit method (|3.5[) . so the proof is complete. □ 
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